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■ Abstract. We will report on the status of a new large scale calculation of thermodynamic quantities 

' in QCD with light up and down quarks corresponding to an almost physical light quark mass value 

and a heavier strange quark mass. These calculations are currently being performed on the QCDOC 
Terafiops computers at BNL. We will present new lattice calculations of the transition temperature 
and various susceptibilities reflecting properties of the chiral transition. All these quantities are of 
1^ . immediate interest for heavy ion phenomenology. 
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£| ; INTRODUCTION AND LATTICE SETUP 

O , 

The lattice action we use is designed for finite temperature Monte Carlo Simulations. 
In the gauge sector we use a (2 x 1)-Symanzik improvement scheme which eliminates 
all cut-off effects of order C(a 2 ), where a is the lattice spacing. For the fermions we 
use the staggered fermion formulation. On top of that we add an improvement term 
which restores the rotational symmetry of the free quark propagator on the lattice up to 
order &{p A ) in the momentum p 111]. In order to improve the flavor symmetry, which is 
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violated in the staggered fermion formulation, we smear the gauge field which is used 
in the standard part of the fermion action. The procedure involves either only three link 
terms (p4fat3) or up to seven link terms J2j](p4fat7). While the p4fat3 action was used 
for thermodynamical calculation earlier JJlj], the combination of the p4-term and seven 
link smearing is used here for the first time in thermodynamic calculations 

We perform simulations with three degenerate quark flavors as well as with 2 light and 
one heavy quark flavor. The lattices have temporal extent N t = 4, 6, which corresponds at 
the critical temperature (T c ) to a lattice spacing of a ~ 0. 13 fm and 0.22 fm respectively. 
The lattice extent in spacial direction is N s = 8,16,32. To determine the scale, we 
perform zero temperature simulations on 16 3 x 32 lattices. These calculations are being 
performed on the QCDOC Terafiops computers at BNL. 



THE CRITICAL TEMPERATURE 

The observables we calculate at present are the gauge action, the Polyakov loop and the 
chiral condensate. For each of these quantities we calculate also its susceptibility. Using 
the multi-histogram reweighting technique [7] we combine the statistics of Monte Carlo 
runs at different couplings /3 . Results for the chiral condensate and chiral susceptibility 
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FIGURE 1. The chiral condensate and (left) and the chiral susceptibility (right) for the p4fat7 action 
with three degenerate flavor and N t = 4. Data points are the individual Monte Carlo runs; solid curves are 
from multi-histogram reweighting. 




FIGURE 2. The critical couplings vs. the bare quark masses for different actions and lattice sizes (left) 
and the critical temperature in units of the string tension (a) as a function of the ratio of pseudo-scalar 
and vector meson masses (right). 



of the p4fat7 action in the case of three degenerate flavors and N t = 4 are shown in 
Fig. [U As expected the chiral transition becomes more pronounced for smaller quark 
masses. The peak heights of the susceptibilities increase and a volume dependence 
becomes visible. We define the critical couplings (/3 C ) as the peak positions of the chiral 
susceptibilities. However, we find that the peak positions of all susceptibilities agree 
within our statistical accuracy. In Fig. Etleft) we show a comparison of /3 C for the two 
actions and temporal lattice sizes analyzed so far. The solid lines are fits to a power law. 
We find that the separation between the N t =4 and N t = 6 critical lines is much smaller 
for the p4fat7 action than for the p4fat3 action. In fact, for small quark masses the critical 
couplings for N t —4 and 6 almost coincide which suggests the presence of a nearby bulk 
transition in the parameter space defined by the p4fat7 couplings. This also leads to the 
onset of first order transitions at pion mass values larger than those found for the p4fat3 
action. All this makes the p4fat7 action less suitable for studies of thermodynamics 1 



We note that we have used a p4fat7 action with tree level coefficients and have not introduced tadpole 
improvement factors as it is done for the asqtad action @|. 



as the extrapolation to the continuum limit will become more difficult, i.e. will require 
larger N t . 

The two (2 + 1) -flavor lines shown in Fig.Etleft) are lines at constant strange quark 
mass. The lines corresponds to am s = 0.065 and am s = 0. 1 for the lower and upper curve 
respectively. A bare strange quark mass of am s = 0.065 approximately corresponds to 
the physical strange quark mass. 

At each critical coupling /3 C we perform a zero temperature scale setting calculation on 
a 16 3 x 32 lattice, where we calculate hadron masses and heavy quark potentials. From 
the potential we extract the Sommer scale tq as well as the string tension a by fitting the 
lattice data with the usual Ansatz V(r) = —a/r + or + const. We show some preliminary 
results of the critical temperature in units of the string tension in Fig. fright). The quark 
mass dependence of the critical temperature for p4fat3, N t = 4 and nf = 3 is consistent 

with the earlier calculations 0]. The chiral extrapolation of T c / \f((j) with (mps/my) 2 
yields results consistent with earlier findings of T c / \/~((j)\m=o = 

0.407(15) 0]. The 

chiral extrapolation for p4fat7 gives smaller values: T c /^/(o)\ m= o < 0.4 and < 0.35 
for N t = 4 and N t = 6 respectively. Furthermore, the N t dependence of the transition 
temperature is much stronger for p4fat7 than for p4fat3. This again indicates a possible 
bulk transition in the vicinity of the chiral limit of the p4fat7 action. 

OUTLOOK 

The still ongoing determination of the critical temperature is part of a detailed study of 
the thermodynamics of (2+l)-flavor QCD. As soon as we have reliable chiral extrapola- 
tions for N t = 4 and N t = 6 we will add a N t =S calculation at selected quark mass values 
and try to perform the extrapolation to the continuum limit. Moreover we will perform 
a study of the equation of state with an almost realistic quark mass spectrum. 
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